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Optical excitation of the sequential super molecule H2P — ZnP — Q induces 
an electron transfer from the free-base porphyrin (H2P) to the quinone (Q) 
via the zinc porphyrin (ZnP). This process is modeled by equations of 
motion for the reduced density matrix which are solved numerically and 
approximately analytically. These two solutions agree very well in a great 
region of parameter space. It is shown that for the majority of solvents the 
electron transfer occurs with the superexchange mechanism. 

I. INTRODUCTION 

The investigation of photoinduced charge transfer is important both for the de- 
scription of natural photosynthesis [] and for the creation of artificial photoenergy- 
converting devices [] . For experimental realizations of such artificial devices porphyrin 
complexes are good candidates []. Of major interest are those complexes with an ad- 
ditional bridging block between donor and acceptor []. 

Electron transfer reactions can occur through different mechanism []: sequential 
transfer (ST) or superexchange (SE). Changing a building block of the complex [ 
or changing the environment [] can modify which mechanism is most significant. 
To clarify which mechanism is present one sequentially varies the energetics of the 
complex []. This is done by radical substituting the porphyrin complexes [] or by 
changing the polarity of the solvent [] . Also the geometry and size of a bridging block 
can be varied and in this way the length of the subsystem through which the electron 
has to be transfered []. 

SE [] occurs due to coherent mixing of the levels [] and plays a role for any detuning 
of the energy levels []. The transfer rate in this channel decreases exponentially with 
increasing length of the bridge []. When incoherent effects such as dissipation and 
dephasing dominate [], the transfer is mainly sequential [], i. e., the levels are occupied 
mainly in sequential order []. An increase in the bridge length induces only a small 
reduction in the transfer rate []. 

In the case of coherent SE the dynamics is mainly Hamiltonian and can be de- 
scribed on the basis of the Schrodinger equation. The physically important results 
can be obtained by perturbation theory [], most successfully by the Marcus theory 
[]. In case of ST the environmental influence has to be taken into account. The 
more natural description of the relaxation process is based on the density matrix 
(DM) formalism []. The master equation that governs the DM evolution as well as 
the appropriate relaxation coefficients can be derived from such basic information as 
system-environment coupling strength and spectral density of the environment []. 

The main physics of the system can be described by a DM equation which accounts 
for relaxation effects phenomenologically []. The master equation is analytically solv- 
able only for the simplest models []. Most investigations are based on the numerical 
solution of this equation []. However, an estimations can be obtained within the 
steady-state approximation []. Here we perform numerical as well as approximate 
analytical calculations. 

II. MODEL 

We investigate the photoinduced electron transfer in supermolecules that consist 
of sequentially connected molecular blocks, namely donor, bridge, and acceptor. The 
donor (D) is not able to transfer its charge directly to the acceptor (A) because of 
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FIG. 1. Schematic 
view of the energy lev- 
els in the H 2 P-ZnP- 
Q complex taken into 
account in calculation. 
The three states in the 
boxes define the charge 
separation which can 
happen either by ST or 
by SE. 



their spatial separation. D and A can exchange their charges only through B (Fig. 1). 
In the present investigation the supermolecule consists of free-base porphyrin(i/ 2 -P) 
as donor, zinc substituted porphyrin(ZnP) as bridge, and benzoquinone as accep- 
tor []. In each of those molecular blocks we consider only two molecular orbitals, 
the LUMO and the HOMO. Each of those orbitals can be occupied by an elec- 
tron (|1)) or not (|0)). This model allows us to describe the neutral nonexcited 
molecule |1) homo\®) lumo and the following three states of the molecule: neutral 
excited molecule \0)homo\1)lumo, positive ion |0) homo I 0)lumo > and negative ion 
\1) homo\1) lumo- Below Roman indices indicate molecular orbitals (m — - HOMO, 
m—1- LUMO), while Greek indices indicate molecular blocks (/i — 1 - donor, /i — 2 
- bridge, \x — 3 - acceptor). Each of the electronic states has its own vibrational 
substructure. However the time of vibrational relaxation [] is two orders of magni- 
tude faster than the characteristic time of the electron transfer []. Because of this 
we assume that only the vibrational ground states play a dominant role in electron 
transfer. 

One can describe the occupation of an orbital by an electron with the appropriate 
creation operator c+ m = | l)^ m (0\^ m as well as its annihilation c^ m = |0)^ m (l| 
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gives the number of electrons in the molecular block fx. 



For the description of charge transfer and other dynamical processes in the system 
we introduce the Hamiltonian 



H = H S + H B + H 



l SB 



(1) 



where H s characterizes the supermolecule, H B the dissipative bath, and H SB the 
interaction between the two. H s , however, includes the static influence of the envi- 
ronment, namely of the solvent dipoles, which gives rise to a reduction of the energy 
levels, 
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{E e i + E ion ) + V , 



(2) 



The energies E^ m are calculated in the independent particle approximation []. e s 

denotes the static dielectric constant of the solvent. E e i = E/X^V ~~ l) e V(4eo?>) 
describes the energy to create an isolated ion. This term depends on the characteristic 
radius r M of the molecular blocks. E ion = XL EiX^V — — l) e V (^ 7re o r ^) includes 
the interaction between the already created ions. It depends on the distance between 
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hopping term V = E^ *V 



The last contribution to the system Hamiltonian is the 
+ - l ? + K - !) 2 )> wh ich includes the 



coherent hopping between each pair of LUMO V 
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as well as 



the corresponding intensities v^ u . The matrix elements of this operator give nonzero 
contribution only if one of the states has a charge separation. Because there is no 
direct connection between donor and acceptor we assume v± 3 = 0. 



As usual the bath is given by harmonic oscillators with creation and anhilation 
operators and a\. The system bath interaction comprises both irradiative and 
radiative transitions. For t <C 1 — 10 ns the latter one can be neglected . The 
irradiative contribution corresponds to energy transfer to the solvent and spreading 
of energy over vibrational modes of the supermolecule 

H SB = E E K ^ v »M + *x)(y+ + v- v ) , (3) 

where K\^ v reflects the interaction strength between bath mode A and quantum 
transition between LUMO levels of molecules \x and v. 

Initially we use the whole density matrix of system and bath for the description 
of the dynamics. After applying the Markov and rotating wave approximations and 
tracing out the bath modes [] we obtain the equation of motion for the reduced density 
matrix (RDM) 

a = -i/h[H s ,a] + E r M4K^i) + l)((V'^] + [^'^]) 

+«Ki,i)([^>, V~\ + aV-])} , (4) 

where the dissipation intensity Y^ u = 7rK' \ v p{w yXvxjv 1 depends on the coupling K^ u 
of the transition /xl ~> ul and on the bath mode of the same frequency. T^ v depends 
also on the density p of bath modes at the transition frequency and on the 

corresponding coherent coupling v^ u between the system states. n{ui) denotes Bose- 
Einstein distribution. 

For simplicity we introduce a superindex i = {jum}, the intensities of the dissi- 
pative transitions = Tijn(-ujij) between each pair of states, as well as the corre- 
sponding dephasing intensities 7^ = 1/2 J2k(dik + d^j)- Taking these simplifications 
into account one gets 

& u = —i/U^iVijCTji - (TijVji) - E dijCTu + E djiO'jj , (5) 

j i j 

&ij = {-iujij - jij)aij - i/hVijio-jj - an) . (6) 

The simplification is that we do not calculate the system parameters, rather we 
extract them from experimental data. 

III. EXTRACTION OF SYSTEM PARAMETERS 

The porphyrin absorption spectra [] consist of high frequency Soret bands and 
low frequency Q bands. In case of ZnP the Q band has two subbands, Q(0, 0) and 
Q(l, 0). In the free-base porphyrin H 2 P the reduction of symmetry induces a splitting 
of each subband into two, namely Q x (0,0), Q y (0,0) and ^(ljO), Q y (l,0). So the 
emission spectra of ZnP and H^P consist of two and four bands, respectively. Each 
of the abovementioned spectra can be represented as a sum of Lorentzians with good 
precision. It is important to note that the spectra of porphyrin complexes contain all 
bands of the isolated porphyrins without essential changes. We use the lowest band 
of each spectrum. The corresponding frequencies and widths are shown in table [l|. 

On the basis of the experimental spectra we determine Ed*ba = 1.82eV and 
Edb*a = 2.03eV (in CHiCli). The authors of Ref. give the energies of two other 
levels, E D + B - A = 2A4eV and E D + BA - = lA2eV. This allows to calculate E DB + A - = 
1.21eV. The hopping intensity v 23 = v — 2.2meV is calculated in Ref. . On the 
other hand Rempel et al. [] estimate the electron coupling of the initially excited and 
charged bridge states v 12 = V = 65meV. We take the intensity of the intermolecular 
conversions T2i, T23 in range 1 — 10 x 10 11 s _1 []. 



Table 1: Low-energy bands of the porphyrin spectra for CH2CI2 as solvent. 





Absorption 


Emission 




Frequency, eV Width, eV 


Frequency, eV Width, eV 


H 2 P 


i/gb = 1.91 70*0 = 0-06 


z^i = 1.73 7^ = 0.05 


ZnP 


v m = 2.13 700 = 0.07 


z/qi = 1.92 701 = 0.05 



The main parameter which controls the electron transfer in a triad is the relative 
energy of the state D + B~ A. This state has a strong coupling to the solvent that 
changes the energy of the state. The values of the energy E D + B -^ calculated in 
the present model are shown in table for some solvents. In table 2 e s denotes the 
static dielectric permittivity, the optic dielectric permittivity, MTHF 2-methil- 
tetrahydrofuran, and CYCLO denotes cyclohexane. The calculated value E e +b-a — 
2.86eV deviates 15% from the data of Ref. . 



IV. RESULTS 

The time evolution of charge transfer within the supermolecule is described by 
Eqs. (H|) and @. At initial time only the donor state is occupied. The calcula- 
tions were performed with two methods, direct numerical integration and analytic 
approximation. 

For the numerical simulation the eigenvalues and -vectors of the system are cal- 
culated and with these the time evolution of the system is known. The simulation of 
the system dynamics with the parameters determined in the previous section shows 
exponential growth of the acceptor population. Such a behavior can be accurately 
fitted to the formula ^(i) = Ps(oo)[l — exp (— k ET t)}, where Uet — 5 x 10 9 s _1 and 
P 3 (oo) ~ 0.95 for CH2CI2 as solvent. The population of the bridge state does not ex- 
ceed 0.005. This shows that the SE mechanism dominates over the ST for the chosen 
set of parameters. In this case the system dynamics can be described by two values: 
the acceptor population at infinite time P^{oo) and the reaction rate k E T that we 
deduce from the dynamics via the following formula k E r = Psi 00 )/ {Io°[^- ~ Ps(t)]dt}- 

The analytical approach is valid for the kinetic limit t ^> l/^ij- In Laplace-space 
we can replace 1/(^^+7^ + 5) by l/(iuij+jij), where s denotes the Laplace variable. 
This allows to simplify Eqs. ([5]) and (&) and we define a new relaxation operator 

j9ji a jj- 111 this expression the transition coefficients g^ 
contain both, dissipative and coherent contributions 

gij = dij + VijVjijij/lh 2 ^ + 7J)]. (7) 

Assuming the bridge population to be zero allows us to find the dynamics of the 
acceptor state in the form -Ps(t) = P3(oo)[l — exp (— k E Tt)}, where the final population 
P3(oo) and the reaction rate k E T are expressed in terms of the coefficients 

, . 923(912 — 932) D , v. 912923 /, x-1 / Q x 

k ET = g23 H ; , Pi(°o) = ■ (kET) ■ (8) 

921 + 923 921 + #23 

Table 2: Energy of the charged bridge state and transfer rates in different solvents. 



Solvent 


75%CH 2 Cl 2 
+2h%CH z CN 


CH2CI2 


MTHF 


CYCLO 


e s 


15.75 


9.08 


6.24 


2.02 


£00 


2.00 


2.01 


2.03 


2.03 


Ed+b-a, eV 


1.89 


2.86 


3.18 


5.30 


kET, 


3.98 x 10 11 


5.01 x 10 9 


7.94 x 10 8 


3.80 x 10 8 
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FIG. 2. The dependence of the reaction rate (upper row) and final population of the 
acceptor state (lower row) on the parameters V = v±2, v = V23, T = T21, 7 = ^3. Solid 
lines correspond to the numerical solution and dashed lines to the analytical solution. The 
circles show the realistic parameter values for CH2CI2 as solvent. 
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FIG. 3. Dependence 
of reaction rate (upper 
picture) and final accep- 
tor population (lower pic- 
ture) on the energy of 
the bridge state E = 
E D + B ~ A . Solid lines cor- 
respond to the numeri- 
cal solution and dashed 
lines to the analytical so- 
lution. 



V. DISCUSSION 

The following question will now be discussed: How does the mechanism and speed 
of the reaction depend on a deviation of the parameters from the determined values? 
Namely which parameters have to be changed in order to change not only the reaction 
rate quantitatively, but the dominant mechanism of reaction and the qualitative 
behavior of dynamics at all. To answer these questions we calculate the system 
dynamics while varying one parameter at a time and keeping the other parameters 
unchanged. The dependencies of transfer rate Uet and final population Pzipo) on 
coherent couplings V = vu, v = V23 and dissipation intensities T = T 2 i, 7 = r 23 are 
shown in Fig. 2. 

In particular, the decrease of the coherent coupling V induces a quadratic decrease 
of the reaction rate k B T until saturation V ~ lO^ps -1 . Then k^T reaches its lower 
bound and does not depend on V anymore. This corresponds to the crossover of the 
reaction mechanism from SE mechanism to ST. But, due to the big energy difference 



between donor and bridge state the efficiency of this ST is extremely low, i. e., P3 ~> 0. 
The considered variation of the coherent coupling can be experimentally performed 
by exchanging building blocks in the supermolecule. 

The most crucial change in the reaction dynamics can be induced by changing 
the energies of the system levels. As discussed above this can be done by altering the 
solvent. Most important is the relative energy of the bridge state \D + B~A). The 
results of the corresponding calculations are presented in Fig. 3. For high energies of 
the bridge state E D + B - A ^> E D * BA the numerical and analytical results do not differ 
from each other. The reaction occurs with the SE mechanism that coincides with 
the conclusion of Ref. . This is the case for the most of solvents (see table 2). The 
smooth decrease of energy induces an increase of the reaction rate up to the maximal 
value near lps -1 . 

While the bridge energy approaches the energy of the donor state the ST mecha- 
nism starts to contribute to the process. As can be seen in table 2 this regime can be 
reached by the use strong polar solvents. The analytical solution does not coincide 
with the numerical one anymore because the used approximations are no more valid 
in this region. In the case E D + B - A < E D + BA - one cannot approximate the dynamics 
of the acceptor population in the form P 3 ~ [1 — exp (—k B Tt)]- A high value of the 
bridge energy ensures the transition of the whole population to the acceptor state 
\D + BA~). In the intermediate case, when the bridge state has the same energy as the 
acceptor state, the final population spreads itself over these two states ^3(00) = 0.5. 
At even lower bridge energies the population gets trapped at the bridge state. 

We performed calculations for the electron transfer in the supermolecular complex 
H 2 P — ZnP — Q within the RDM formalism. The resulting analytical and numerical 
reaction rates are in good agreement with each other and in qualitative correspon- 
dence with experimental data []. The SE mechanism of electron transfer dominates 
over the sequential one. The qualitative character of the transfer reaction is stable 
with respect to a small variation of the system parameter. The crossover between the 
reaction mechanisms can be forced by lowering the bridge state energy to the energy 
of the donor state. 
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